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A recent study [Tassel et al., Phys. Rev. Lett. 105, 167205 (2010)] has proposed a remarkable 
spin model for (CuCl)LaNb207, in which dimers are ferromagnetically coupled to each other on 
the distorted Shastry-Sutherland lattice. In this model, the intra-dimer exchange coupling J > 
is antiferromagnetic, while the inter-dimer exchange couplings are ferromagnetic and take different 
values, Jx,Jy < 0, in the two bond directions. Anticipating that the highly frustrated character of 
this model may lead to a wide range of behaviors in (CuCl)LaNb207 and related compounds, we 
theoretically investigate the ground state phase diagram of this model in detail using the following 
three approaches: a strong-coupling expansion for small Ja, and Jy, exact diagonalization for finite 
clusters, and a Schwinger boson mean field theory. When \ Jx\,\Jy\ ^ J, the system stays in a dimer 
singlet phase with a finite spin gap. This state is adiabatically connected to the decoupled-dimer 
limit Jx = Jy = 0. We show that the magnetization process of this phase depends crucially on 
the spatial anisotropy of the inter-dimer couplings. The magnetization shows a jump or a smooth 
increase for weak and strong anisotropy, respectively, after the spin gap closes at a certain magnetic 
field. When \Jx\ or \Jy\ > J, quantum phase transitions to various magnetically ordered phases 
(ferromagnetic, coUinear stripe, and spiral) occur. The Schwinger boson analysis demonstrates that 
quantum fluctuations split the classical degeneracy of different spiral ground states. Implications for 
(CuCl)LaNb207 and related compounds are discussed in light of our theoretical results and existing 
experimental data. 

PACS numbers: 75.10.Jm, 75.10.Kt, 75.40.Mg 



I. INTRODUCTION 

The search for exotic quantum states in tvifo- 
dimensional spin systems withJ[i|Ustrated interactions has 
been of great recent interest Jj'EI Among various materi- 
als studied recently, S — ^ layered copper oxyhalidcs 
(CuX)A„_iB„03„+i offer an interesting family of frus- 
trated magnets with rich variety of behaviors. In this 
family, each magnetic CuX layer (with X=Cl,Br) is sand- 
wiched by nonmagnetic layers, forming an ideal two- 
dimensional structure. In each layer, the Cu^+ ions 
form a square lattice of 5* = | spins, and a competi- 
tion between ferromagnetic and antiferromagnetic inter- 
actions is anticipated from the small Curie- Weiss tem- 
peratures (relative to other characteristic energy scales) 
commonly observed in this family. Extensive experi- 
mental investigations have uncovered a collective singlet 
ground state with a spin gap in (CuCl)LaNb207 (Refs. ^- 
g), a collinear stripe magnetic order in (CuBr)LaNb207 
(Ref. and a magnetization plateau at 1/3 of the sat- 
urated moment in (CuBr)Sr2Nb30io (Ref. ||). Further- 
more, chemical substitution was used to observe quantum 
phase transitioris, between magnetic and non-magnetic 
ground states.tHlj It would be interesting to determine 
what kind of spin models can capture the wide variety of 
physics in this family. A frustrated J1-J2 model on the 
square lattice with ferromagnetic Ji and antiferromag- 
netic J2, as initially postulated for this family, does not 
seem to exhibit a spin gap or a magnetization plateau.tHI 

A recent study by Tassel et a/.Oj has proposed a re- 



markable microscopic structure in (CuCl)LaNb207. The 
x-ray and neutron diffraction studies have identified a 
considerable distortion of the Cu-Cl bonds; consequently, 
the unit cell is doubled along a and b axes, as shown in 
Fig. 0(a). Based on the obtained crystal structure, a den- 
sity functional calculation was carried out to construct 
the microscopic spin model. The dominant interaction 
was then found to be the fourth-neighbor antiferromag- 
netic exchange coupling J, which pairs spins into dimer 
singlets shown as thick lines in Fig. |^(a). The calculation 
also suggested that the next leading couplings are ferro- 
magnetic couplings (denoted by and Jy in this paper) 
between these dimers, which are shown by solid and bro- 
ken lines in Fig. ^ (a) . Remarkably, the resultant interac- 
tion network has the structure of the distoxtcd Shastry- 
Sutherland latticco as shown in Fig. |l|(b).li3 

In this paper, we study a spin-i model of coupled 
dimers on the distorted Shastry-Sutherland lattice, as 
the simplest starting point for understanding the inter- 
play of competing interactions in (CuCl)LaNb207 and 
related compounds. We consider the Heisenberg model 
in a magnetic field. 



(1) 



where runs over all the bonds in Fig. 0(b), and 

Jij = J, Jx, and Jy for diagonal, horizontal and verti- 
cal bonds, respectively. We are primarily concerned with 
the case of antiferromagnetic J > and ferromagnetic 
Jxt Jy < 0; however, results for other signs of Jx and Jy 
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FIG. 1: (Color online) (a) Schematic picture of a single mag- 
netic layer of (CuCl)LaNb207. Thick solid (J), thin solid 
(Jx), and broken (Jy) lines indicate the three major couplings 
revealed in the electronic structure calculation of Ref. uA. 
Other interactions J2b an d Jl shown by dotted lines will 
also be considered in Sec. VII. (b) Topologically equivalent 
Shastry-Sutherland picture of the three major couplings J, 
Jx, and Jy. A dotted green plaquette indicates a unit cell. 
The four sublattices are labeled as A, B, C, and D. 



are presented alongside for comparison. A common view- 
point for dealing with coupled dimer systems is to regard 
the triplet excitation at each dimer as a particle (called 
a "triplon") andJ-o-describe the system as a Bose gas 
of such particles. E£lEZl A notable feature of the Shastry- 
Sutherland lattice is the strpng^uppression of the triplon 
hopping due to frustration.E§E3 In the antiferromagnetic 
case with J > and Jx = Jy > 0, the localized nature 
of the triplons gives rise tcussaxious fractional plateaux in 
the magnetization process,c2rE3 which are experimentally 
observedE^ in SrCu2 (603)2. Frustration also exist in the 
case of our interest, with ferromagnetic Jx, Jy < 0. The 
consequence of frustration in the ferromagnetic case has 
not been addressed in previous studies. 

Our main results for the ferromagnetic case Jx, Jy < 
are as follows. When | |</y| ^ J, the system stays in 
a dimer singlet phase with a finite spin gap. This state 



is adiabatically connected to the decoupled dimer limit 
Jx = Jy = 0. The magnetization process of this phase 
depends crucially on the spatial anisotropy of the inter- 
dimer couplings: the magnetization shows a jump and a 
smooth increase for weak and strong anisotropy, respec- 
tively, after the spin gap closes at a certain magnetic field. 
When \Jx\ or \Jy\ > J, the spin gap of the dimer sin- 
glet phase closes, and quantum phase transitions to vari- 
ous magnetically ordered phases (ferromagnetic, coUinear 
stripe, and spiral) occur. It is demonstrated that quan- 
tum fluctuations split the classical degeneracy of different 
spiral ground states. The detailed phase diagrams of the 
classical and quantum models are constructed. These re- 
sults are based on a strong-coupling expansion for weak 
Jx and Jy, exact diagonalization for small clusters, and 
a Schwinger boson mean field theory. Using the existing 
data on the magnetization process and the triplon band- 
width in (CuCl)LaNb207, and comparing with the cor- 
responding theoretical results, we provide a consistency 
check for the appropriateness of the proposed spin model 
(I). 

The rest of the paper is organized as follows. In 
Sec. we present the classical and quantum phase di- 
agrams of the model (|^) and summarize the main re- 
sults of the paper. In Sec. |ll|, we discuss the details of 
the computations for the classical model and the result- 
ing ground states. In Sec. IV, we perform the strong- 



coupling expansion of the quantum model to derive an 
effective low-energy Hamiltonian. This effective model 
will be very useful in determining the global phase dia- 
gram of the quantum model and in estimating relevant 
physical quantities. In Sec. the numerical analysis 
of the quantum model based on the exact diagonaliza- 
tion is presented. Some parts of the phase boundaries 
are determined quite accurately within this approach. In 
Sec. VI, the Schwinger boson mean field theory analy- 
sis of the model is reported. These three techniques are 
collectively used to understand the quantum phase di- 
agram. In Sec. VII, we compare these theoretical re- 



sults on the quantum model with experimental da ta on 
(CuCl)LaNb207. We conclude the paper in Sec. VIII. 
Appendix ^ explains an exact solution available in the 
isotropic case Jx — Jy In Appendix ^ the Hartree vari- 
ational approach is discussed. In Appendix the sat- 
uration field in the magnetization process is determined 
exactly by considering a single-magnon excitation from 
the fully polarized state. 



II. PHASE DIAGRAMS AND MAIN RESULTS 

In this section we present the classical and quantum 
phase diagrams of the model (Q) in Figs. |[ ||, and ||, and 
summarize the main results of the paper. We assume J > 
and h> throughout the paper. With the assumption 
of J > 0, frustration at the classical level occurs only 
when JxJy > 0. When we discuss the isotropic case 
Jx — Jy, we denote these two parameters by J'. 
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FIG. 2: (Color online) Classical phase diagram of the dis- 
torted Shastry-Sutherland model (^) with h = 0. Green 
dashed lines at Jx = 0, Jy — 0, Jx = Jy are guides for eyes. 



states are degenerate, apart from the trivial degeneracy 
associated with the global 0(3) rotation of spins. One 
ground state is the "a;-spiral" in Fig. |6|(c), where spins 
rotate by a uniform angle Qx in the x direction and al- 
ternating angles ±Qy in the y direction. The other is the 
"y-spiral" , which is defined by interchanging the roles of 
X and y directions. Remarkably, in the anisotropic case 
Jx 7^ Jy, the X- and y-spirals are not related to each other 
by any symmetry operation. Their degeneracy therefore 
comes from the particular geometry of the lattice. 

We note that the classical phase boundaries in Fig. ^ 
are symmetric with respect to the sign flips Jx,y — > —Jx,y 
This is because the signs of Jx and Jy can both be 
flipped by reversing the spins on the A and D sites of 
all units cells, as seen in Fig.0(b). This transformation 
is not allowed in the quantum case, since the simultane- 
ous reversal of all the spin components (S^ , , S^) ^• 
{—S^,—S'^,—S^) changes the commutation relations of 
spin operators. The quantum phase diagram presented 
next therefore depends on the signs of Jx and Jy. 



B. Quantum phase diagram 



A. Classical phase diagram 

The classical spin model is defined by replacing all 
the spin operators in Ec^. (P by classical 0(3) vectors 
of length S. Assuming h — for simplicity, the ground 
state phase diagram in Fig. || is determined exactly for 
different signs of Jx and Jy. The details of our analysis 



will be presented in Sec. III. Below we summarize the 



main characteristics of each phase. 

Collinear stripe phases. These phases appear natu- 
rally in the unfrustrated regions JxJy<0 and penetrate 
into some parts of the frustrated regions JxJy > 0. The 
spin configuration has a propagation vector of {kx, ky) := 
(k ■ x,k ■ y) = (0, tt) or (tt, 0), where x and y are defined 
in Fig. |l|(b). In the ground state with (0,7r), for exam- 
ple, the system forms stripes of up or down spins running 
along the x direction, while up and down alternate in the 
y direction. 

Neel phase. This phase appears for large positive 
Jx/J and Jyl J , where there is a minimal effect from 
the J coupling. The spin configuration is collinear and 
has a propagation vector of {kx, ky) = (tt, tt), as usual for 
the square lattice antiferromagnet. 

Ferromagnetic (FM) phase. This phase appears for 
large negative Jx/J and Jy/J, where ferromagnetic Jx 
and Jy dominate over antiferromagnetic J. The spins arc 
all aligned in the same direction and have a propagation 
vector of {kx,ky) ~ (0,0). 

Spiral phases. These phases appear when the mag- 
nitudes of J, Jx, and Jy are comparable. This is the 
case where the effect of the frustration is most promi- 
nent. The ground state is a coplanar spiral state with an 
incommensurate propagation vector, as originally found 
in the isotropic case Jx = Jy in Ref. A notable fea- 
ture of the spiral phase is that two kinds of spiral ground 



has been stud- 



Thc quantum model with S = 
led intensively in the isotropic antiferromagnetic case 
Jx = Jy{= J')- A remarkable feature of this model is 
that for /i = and < J'/ J < 1/2, the ground state 
of the Hipdel is exactly given by the product of dimer 
singletalj 



«'ds> = 



n 

R 



s)r, 



(2) 



where R labels a dimer and |s) r is the singlet state on the 
dimer. This solution can be extended to the ferromag- 
netic region — 1 < J'/ J < 0, as described in Appendix ^. 
A notable difference between the ferromagnetic- and 
antiferromagnetic- J' cases occurs when a magnetic field 
h is applied. In the antiferromagnetic- J' case, the pafiftn 
nctization process shows various fractional platcaux,E3L3 
which can be viewed as density wave formations of triplon 
excitations. In the ferromagnetic- J' case, the magnetiza- 
tion process shows a jump from the dimer singlet state 
(^ to the fully satu rated state, as in Fig. ^(a). As will 
be discussed in Sec. [V C, this can be viewed as a conse- 



quence of the phase separation of triplons. 

Our major interest is in the ferromagnetic case of 
Jx, Jy < with spatially anisotropic couplings. For com- 
parison, we will also include some results for other signs 
of Jx and Jy. The ground state phase diagrams are pre- 
sented in Figs. I and |. Figure 
coupling expansion for weak Jx and Jy (Sec 
exact diagonalization for finite clusters (Sec. [Vj^ 
I'^a-U'^yl ^ J, the system stays in a dimer singlet phase 
with a finite spin gap, which is adiabatically connected 
to the decouplcd-dimer limit Jx ^ Jy ^ Q. This phase 
is split into three regions (I, II, and III), in terms of the 
behaviors of the magnetization processes, as shown in 



is based on a strong- 
]) and 
When 
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FIG. 3: (Color online) Phase diagram of the distorted 
Shastry-Sutherland model (0) in the spi n- ^ case, determined 
by the strong coupling expansion (Sec. fvp and exact diago- 
nalization (Sec. ^) analyses. Green dashed lines at = 0, 
Jy = 0, Jx = Jy are guides for the eye. The dimer sin- 
glet (DS) phase is divided into three regions, I, II, and III, 
which are characterized by the magnetization processes in 
Fig. ^. Red broken lines around the origin indicate the re- 
gion boundaries determined from the the first-order effective 
Hamiltonian (see Fig. ^). The square and circular symbols 
are based on exact diagonalization for the number of spins, 
Ns = 16, 20, 6x4. The three different symbol sizes are in the 
order of Ns- The classical ferromagnetic phase boundary is 
superposed on the square symbols, showing good agreement. 
The diamond symbols are the boundaries between the DS- 
II regions and the stripe phases, determined by the fidelity 
susceptibility analysis in Fig. |l^. Narrow spiral phases may 
appear between the DS-I region and the stripe phases for large 
\Jx\/J or \ Jy\/J. For this reason, the DS-stripe phase bound- 
aries (diamond symbols) are not calculated beyond Jx/J or 

Jy/J^ -1.3. 
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FIG. 4: Sketches of magnetization processes in the dimer 
singlet phase: (a) a jump, (b) a smooth increase, and (c) a 
plateau at 1/2 of the saturated moment Ms. 
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FIG. 5: Phase diagram of the distorted Shastry-Sutherland 
model (|l|), determined by a Schwinger boson mean field the- 
ory (Sec. ^) for S — 0.5. Solid lines indicate second- 
order transitions, while broken lines indicate first-order tran- 
sitions. The isolated dimer singlet (DS) state is found for 
small |Jii|,|Ji/| < 0.5 J, while the ferromagnetic state stabi- 
lizes for larger \Jx\ and \Jy\. The stripe state expands from 
the classical case to fill \Jx\ or \Jy \ < 0.55 J. The spiral state 
exists in the middle. The a;-spiral appearing for \Jx\ < \Jy\ 
has an incommensurate long-range order in the x direction as 
in Fig. ^. The y-spiral similarly appears for \Jx\ > \Jy\. The 
dotted line separates the two spiral phases. 



Fig. ^. In the DS-I region, the magnetization M shows a 
jump as in Fig. ^(a). In the DS-II region, the magnetiza- 
tion M smoothly increases towards saturation after the 
spin gap closes at a certain magnetic field, as in Fig. ^(b). 
Finally, in the DS-111 region, the magnetization process 
shows a plateau at 1/2 of the saturated moment Ms Jr^W^ 
more plateaux can appear near the isotropic caseE2rE2l 
Jx = Jy)- When \Jx\/J or \Jy\/J > 1, the spin gap of 
the dimer singlet phase closes, leading to various mag- 
netically ordered states. The characteristics of the ferro- 
magnetic and coUinear stripe phases are rather similar to 
those in the classical model. The transition lines from the 
dimer singlet phase to these two phases are determined 
by exact diagonalization with relatively good accuracy 
(square and diamond symbols in Fig. 

The exact diagonalization for small systems, however, 
is not adequate to find incommensurate spiral phases that 
may intervene between the stripe, dimer singlet, and fer- 
romagnetic phases. In order to get insight into this ques- 
tion, we performed a Schwinger boson mean-field analysis 
(Sec. VI). In this approach, the spin magnitude S can be 
varied freely, and one can discuss how the classical phase 
diagram in the limit S* — > oo (Fig. |^) changes as quantum 
fluctuations are gradually taken into account. The phase 
diagram for S = 0.5 is shown in Fig. 0. In this result, 
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wc find that the spiral phases do appear in some parts of 
the phase diagram, while the presence of various magnet- 
ically ordered phases are somewhat exaggerated, as ex- 
pected in mean-field methods. Providing an interesting 
difference from the classical case, quantum fluctuations 
split the classical degeneracy of the two spiral states, fa- 
voring the X- or y-spirals for | J^; | < | Jj, | and \Jx\ > | Jy | , 
respectively. Although the mean-field method is not use- 
ful in discussing the precise locations of the phase bound- 
aries, it is natural to expect that these spirals appear in 
narrow regions between the DS-I region and the stripe 
phases in Fig. |[ particularly for large \ Jx\/J or \ Jy\/ J. 




C. Relation to the experiments on (CuCl)LaNb2 07 

In experiments on (CuCl)LaNb207; the magnetization 
increases smoothly beyond a critical magnetic field un- 
til it reaches the saturated value.cl This suggests that, at 
least within the J-Jx-Jy model, the ground state should 
belong to the DS-II region as described in Fig. || and 
there should be substantial anisotropy in Jx and Jy . Con- 
sidering the magnetization datan and th|C-|measured en- 
ergy range Ae of the triplon excitations,li3 and compar- 
ing them with theoretical results, we find that within the 
J-Jx-Jy model, one of Jx and Jy may be ferromagnetic 
while the other may be antiferromagnetic. However, de- 
pending on how to interpret the existing experimental 
data, the energy range of the triplon excitations can be 
larger; if this is the case, Jx and Jy can be both ferro- 
magnetic. We plot the triplon dispersions for the cases of 
(a) ferromagnetic Jx{< 0) and antiferromagnetic Jy{> 0) 
and (b) both ferromagnetic Jx,Jyi< 0) (Fig. 14). These 



plots can be used to test the J-Jx-Jy model further and 
to determine the signs of Jx and Jy , when more detailed 
information of the triplon excitations is provided from ex- 
periments. We also discuss the effects of other exchange 
couplings J26 and J4 shown in Fig. |l|(a), which are also 
contained in the model of Rcf. 13, These couplings induce 



oscillating behaviors in the triplon dispersions, which can 
be used as fingerprints of their existence. The details of 
these analyses are presented in Sec. VIL 



III. CLASSICAL GROUND STATE 

In this section we describe the exact solution for the 
classical ground state. The phase diagram is prese nted 
in Fig. ^, with the phases summarized in Sec. II A. For 
simplicity, we restrict ourselves to the case of /i = and 
J, Jx, Jy > 0. The case of J > and Jx, Jy < can be 
treated in parallel by applying the spin reversal trans- 
formation on the A an d D sites, as explained in the last 
paragraph of Sec. [LA. In the case of JxJy < 0, the sys- 
tem is not frustrated and the ground states are naturally 
determined as collinear stripe states. 

To find the ground state of the classical model with 
J, Jx, Jy > 0, wc decompose the Hamiltonian into trian- 



FIG. 6: (Color online) (a) A triangular part of the Hamilto- 
nian. (b) Sketch of how to solve the three-spin problem. A 
coplanar ground state is obtained by forming a triangle using 
three spins, (c) A classical spiral ground state with a pitch 
angle Qx in the x direction ("s-spiral"). The number at each 
site indicates the angle (in the xy plane) of the spin, where 
Qx and Qy axe given by Eq. 



gular parts and solve a single-triangle problem shown in 
Fig. I (a). Here the diagonal coupling J is divided by 
two, since it is shared by two neighboring triangles. By 
rewriting the three-spin energy E/\ as 



„ JxJyJ 
= -. 



-J- + —r + ^] + const., (3) 

Jy J Jx / 



we find that the ground state of E/^ is obtained by mini- 
mizing the length of the vector Si/Jy + 23-2 /J + Sy,/ Jx- 
When the vector lengths 2S/J, S/Jx, and S/Jy have 
comparable magnitudes, the three vectors can form a 
triangle so that the above sum vanishes, as shown in 
Fig. 0(b). This leads to a coplanar spin configuration 
as shown in Fig. ^(a), where the spins rotate counter- 
clockwise by Qx and then by Qy when moving between 
sites 1 — >■ 2 — >■ 3. Another configuration where the spins 
instead rotate clockwise by the same angles also gives a 
ground state. Here Qx and Qy satisfy the relations 



J, 



which are solved as 



^ smQx = ^ sinQy, 

Jy <J x 

\- COS Qx + -f- COS Qy = - 



Jx 
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Using a ground state of locally on every triangle, one 
can construct a ground state of the whole lattice as in 
Fig. ||(c). In this state, the spins rotate by a uniform 
angle Qx in the x direction, and by alternating angles 



±Qy in the y direction. We call this state the "x-spiral." 
Similarly, one can construct the "y-spiral" , where the 
spins rotate by a uniform angle Qy in the y direction 
and by alternating angles ±Qx in the x direction. As 
explained in Sec. O, the degeneracy of these states does 
not result from symmetry, and thus it should be regarded 
as a consequence of the particular geometry of the lattice. 

When 2/ J > 1/ J^^ (with J, J^, Jy > 0), the three 
spins of i?A can no longer form a triangle as in Fig. ^(b) . 
Instead, Si and S3 align antiparallel to S^2, forming a 
coUincar configuration with Qx = Qy — tt. This leads to 
a Neel ground state of the whole lattice. Similarly, for 
1/Jx > '2,/J+ l/Jy and 1/ Jy > 2/ J + 1/Jx, one obtains 
coUinear stripe ground states with propagation vectors 
(0,7r) and (7r,0), respectively. 

The phase diagram obtained above for the antiferro- 
magnetic quadrant Jx, Jy > can be mapped onto the 
ferromagnetic quadrant Jx, Jy < by applying the spin 
reversal mentioned in the last paragraph of Sec. II A. The 



Neel ground state maps onto the ferromagnetic state. 
The (0,7r) stripe state maps onto the (tt, 0) stripe state. 
A spiral state with Qx, Qy > tt/2 maps to a spiral state 
with Qx, Qy < 7r/2. These arguments complete the phase 
diagram for all possible signs of Jx and Jy, as shown in 
Fig. I 
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FIG. 7: Examples of virtual processes in the second-order 
perturbation theory, which lead to correlated hopping terms 
(a) (1 — Mi?)&Jj_|_j,]^foH-i-e2 and (b) nii6^_|_gj^feH+e2, and a three- 
body interaction term (c) (1 — n_R)n_R+ein_R+e2. In each pro- 
cess. Hi acts first on R and i?-|-ei, as indicated by the dotted 
lines, to yield the excited state. It then acts on R and R-\- e2 
to yield the final state of the process. 



IV. STRONG-COUPLING EXPANSION 



In this section we analyze the spin-i 
means of a strong coupling expansion.E2ln3 

J y 



el (0) by 
In this 
0, where 



approach, we start from the limit Jx 
dimers are decoupled from each other. We then pertur- 
batively include the effects of Jx and Jy and deri ve an 
effective Hamiltonian of the model (j^). In Sec. fVA. 



we obtain the effective Hami ltonia n u p to s econd order 



in Jx and Jy. Then, in Sees. [VB and |lV C , we analyze 



the first- and second-order effective Hamiltonians, respec- 
tively, to deduce the physical properties of the original 
Hamiltonian H. The first-order Hamiltonian is exactly 
equivalent to the XXZ model on the square lattice. Us- 
ing the known results on the XXZ model, we map out a 
qualitative phase diagram of the model (|]). The second- 
order Hamiltonian provides more accurate estimations of 
physical quantities and phase boundaries than the first- 
order one. The obtained phase boundaries show a re- 
markable agreement with the exact diagonalization re- 
sult of Sec. 0. The s pin gap and the triplon band width 
derived in Sec. IV C| will be use d to fit the experimental 
data of (CuCl)LaNb207 in Sec. |vh1 



Effective Hamiltonian 



We start from the limit Jx = Jy = 0, where dimers 
are decoupled from each other. We derive an effective 
Hamiltonian of the model (|l|) by perturbatively includ- 
ing the effects of Jx and Jy. We label each dimcr by its 



center position R as in Fig. |l|(b); the two spins on the 
dimer are then denoted by Sri and 5'h2- When h = 0, 
the eigenstates of each dimer consist of a singlet |s) and 
a triplet where /x = 0, ±1 represents the z compo- 

nent of the total spin, 5*^^ -I- 5*^2 ■ When a magnetic field 
/i > is applied, the degeneracy of the triplet is split, 
with = Itt) having the lowest energy of the three. 
Then it is useful to focus on the low-energy sector of the 
Hilbcrt space which consists of \s) and \t+i). We ap- 
ply degenerate perturbation theory to derive an effective 
Hamiltonian in this restricted Hilbcrt space. Specifically, 
we adjust the magnetic field to h = J so that the singlet 
|s) and the up-polarized state become degenerate. 

The unperturbed Hamiltonian is 



H, = J2 [JSm ■ Sr2 - J {Si,, + S^R^)] . 



(6) 



R 



The remaining part of the Hamiltonian is 
Hi^Y. [•^-^■^ ■ ^^+- + •^vS^ ■ S,+y -{h- J)S:] , (7) 

i 

where i runs over all the sites. We treat Hi as a per- 
turbation, assuming that the coefficient of each term is 
sufficiently smaller than J. Note that the coefficient of 
the Zecman term in Hi is /i — J (not h) since the mag- 
netic field of magnitude J is included in Hq] therefore, 
the perturbation theory is most accurate around the in- 
termediate magnetic field h = J. The ground state of Hq 
is macroscopically degenerate with a degeneracy of 2^^""^ , 
where N^^ is the number of unit cells in the system. Let 
Vq be the subspace of the Hilbcrt space spanned by these 



7 



states. Regarding and \s) on each dimcr R as the 

presence and the vacancy of a particle ("triplon"), we 
treat the system as a hard-core boson gas on the square 
lattice of dimer centers. On each dimer R, we define a 
creation operator 6^ = |t+i)(s| and a number operator 
nR = Degenerate perturbation theory up to 

second order yields an effective Hamiltonian 



= Eo 



HI 



off 



h: 



eff 



(8) 



where Ei^ 



}N^cJ is the ground state energy of Hq, 



and H^^ and iJj*^ are the first- and second-order contri- 
butions, respectively. 

The first-order contribution H^^ is calculated by pro- 
jecting Hi onto the ground-state manifold Vq of Hq: 



Hf = 

{RR') 



^{bRbR' + h.c.) + '^nRiiR, 



(9) 



R 



where {RR!) runs over all the nearest- neighbor pairs of 
dimers, and 



Jx ^ Jy 



(10) 



The Jx and Jy couplings produce the triplon hopping 
and interaction terms between nearest-neighbor dimers. 
The two couplings contribute with opposite signs to the 
hopping term, so the coefficient is proportional to J_. 
This leads to a suppression of the boson hopping near 
the isotropic case Jx = Jy, which is a—jpjnarkable fea- 
ture of the Shastry-Sutherland lattice.taEiJ The interac- 
tion term is repulsive or attractive depending on the sign 
of J+. The Zeeman term /i — J in iJi plays the role of 
the chemical potential of the hard-core boson gas. Since 
this term is diagonal in the particle number basis {tir}, 
it contributes only to the first-order Hamiltonian H^^ . 



To calculate the second-order contribution H2^ , we 
take into account various processes where a state in Vq is 
virtually promoted to an excited state of Hq by the opera- 
tion of Hi , and then comes back to Vq by Hi again. Such 
virtual processes give rise to correlated hopping terms 
as well as three-body interaction terms, as illustrated in 
Fig. 1^. The resulting form of is 



H 



off 



2J ■ 

{RR') 



Jl Jl\^ 



j+j- 

4J 



4 

X! X! ~ nR)nR+e,nR+e^^^ 4- riR+e^ {b^j^bR^e^ + h.c.) 



R. u=l 

'Jl 



u J2 



ReVi 



j2 j2 \ ^ f J J J'^ \ 

4j"^ - - '^«) ) i^R+e2bR-e2 + h-C) + J] ( (-l)'^^^"^ " " ) (^Ue„ ^ 

J2 j2 
- -^{nR+el - nR-el){b\fbR+e2 ~ &iJ&_R-e2 + h.C.) - -^{nR+e2 " 'ra_R-e2 ) (&lj&fl+el " ^j^^H-el + h.C. 

[ei ^ 62, J- ^ - J-] . 



_i +b.c.) 



(11) 



Here, the vectors ei and 62 are defined as in Fig. and B. First order: qualitative phase diagram 

63,4,5 are defined via the relations ei = —63 = 65 and 
e2 ~ —e^. Vi and 2?2 refer to the sets of dimers oriented 
in the directions of ei and 62, respectively. By setting 
Jx = Jy, the effective Hamiltonian presented hcre|-cpin- 
cides with the one derived by Momoi and Totsuka.El 

By identifying the hard-core boson operators with the 
spin-i operators via — and ur = s|j + 1/2, the 
first-order effective Hamiltonian Hl^ is equivalent to the 
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XXZ model on the square lattice, 



{RR') 



R 



with 



K. 



J+ 



xy — ^J-: Kz — h ~ h — J — Jj^. 



(12) 



(13) 



Note that this model has a higher symmetry than the 
original model; the unit cell has been reduced from two 
dimers to a single dimer. The phase diagram of —this, 
model has been studied in detail in the literature .EZro 
By changing h from large negative to large positive val- 
ues, the ground state of Hf^ changes between the fully 
down-polarized {{sr) = —1/2) and fully up-polarized 
{{sr) ~ +1/2) states. In the original model (0), this 
corresponds to a change between the dimer singlet state 
(^ with magnetization M — and the fully polarized 
state with M = Ms- The process of this change can be 
smooth or sudden; the details depend on the value of the 
XXZ anisotropy Kz/\Kj;y \ of Hf^ , and are classified into 
three cases (I)- (III), described below. We note that the 
sign of Kxy is not crucial in the present argument, since 
it can be flipped by applying a gauge transformation to 
or by exchanging the roles of Jx and Jy in H. 

(I) Ferromagnetic Ising case Kz/\Kxy\ < —1. At 
h = 0, the fully up- and down-polarized states are degen- 
erate in the ground state of Hf^ . Therefore, the magne- 
tization of the effective XXZ model shows a jump from 
(sr) ~ —1/2 to +1/2 as h changes its sign from nega- 
tive to positive. This corresponds, in the original model, 
to a jump of the magnetization M at an intermediate 
magnetic field he = J + J+ as in Fig. ^a). 

(II) XY case -1 < Kz/\Kxy\ < 1. The ground state 
of Hf^ at /i = is given by a ferromagnetic or anti- 
ferromagnetic state in the xy plane for K^y < and 
Kxy > 0, respectively. This XY ordered state changes 
to a canted state by a magnetic field h and the 
magnetization (sr) varies smoothly between —1/2 and 
+1/2 in the range ~2{\Kxy\ +Kz)<h< 2{\Kxy\ + Kz). 
Therefore, in the original model, the magnetization M in- 
creases smoothly from zero to the saturation in the field 
range between hd = J— 2|J_| and hc2 = J-I-2(J_|_-|-| J_|), 
as in Fig. ^(b). The canted ferromagnetic and antiferro- 
magnetic orders in the effective model give rise to trans- 
verse magnetizations in the original model. Specifically, 
in the field range hd < h < hc2, the transverse compo- 
nent shows a stripe order with a propagation vector of 
{kx,ky) = (0,7r) and (7r,0), for < Jy and Jx > Jy, re- 
spectively. The connection between the magnetic orders 
of the effective and original models is understood most 
easily from a variationalground state approach, which is 
presented in Appendix 

(III) Antiferromagnetic Ising case Kz/\Kxy\ > 1. In 
contrast to the XY case above, the ground state of Hf^ 
at /i = is given by an antiferromagnetic state polarized 
along the z axis, featuring gapped excitations. When h is 



Stripe 




1 1-5 J /J 



FIG. 8: Qualitative phase diagram of the spin-i model 
derived from the first-order effective Hamiltonian H^^ of the 
strong-coupling expansion. In the DS-I and DS-II regions, 
the magnetization shows a jump and a smooth increase as 
in Figs, ^(a) and (b), respectively. In the DS-III region, at 
least one magnetization plateau is expected to appear at 1/2 



of the saturated moment as in Fig. 
be expected near the isotropic case. 



More plateaux may 



increased above zero, the magnetization (sr) stays zero 
up to certain critical h > and then shows a jump in the 
transition to the canted antifcrceunagnetic state (known 
as the spin-flopping process). EZrEa A similar thing occurs 
when decreasing h below zero. Therefore, in the original 
model, the magnetization process shows a plateau at 1/2 
of the saturation, accompanied by jumps at the edges, as 
in Fig. ^(c). The expressions of hd and hc2 are the same 
as the XY case (II) above. 

The present argument allows us to map out a qual- 
itative phase diagram as in Fig. |^. The dimer singlet 
phase appearing for | J^, | , | Jy | < J are divided into three 
regions (I, II, and III) in terms of the behaviors of the 
magnetization processes shown in Fig. ^. The boundaries 
between the regions are given by Kz = iKxy from the 
above argument. In principle, the effective Hamiltonian 
can give reliable results for H only when Jx and Jy 
arc sufficiently smaller than J. Beyond this region, we 
nonetheless use Hf^ to investigate possible instabilities 
of the dimer singlet state. By increasing \Jx\ and \Jy \ in 
the DS-I region, the critical field he = J + J+ goes to 
zero, indicating a transition to the ferromagnetic phase. 
Similarly, by increasing \Jx\ or \Jy \ in the DS-II regions, 
hd — J — 2|J_| goes to zero, leading to the collinear 
stripe phases. The boundaries between the ferromag- 
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netic and stripe phases are determined by the condition 
hc2 = J + '2{J+ + I J_|) = 0. The obtained phase diagram 
in Fig. H quahtatively agrees with Fig. || obtained by ex- 
act diagonahzation. However, some aspects of the phase 
diagram are not captured in this approach. For example, 
the spiral and Neel phases (expected from the classical 
analysis of Sec. [II and the Schwingcr boson analyses of 
Refs. n 



30y3l| and Sec. VI) do not appear in this approach. 
These phases are beyond the scope of the present pcrtur- 
bative analysis. 

In passing, we note that the phase diagram in Fig. ^ 
can also be obtained using a Hartree variational state 
approach, which is presented in Appendix H. 



C. Second order: estimations of physical quantities 

We now move to the analysis of the second-order ef- 
fective Hamiltonian . On one hand, this allows us to 
make more accurate estimations of physical quantities, 
such as the spin gap, compared to the first-order case 
above. On the other hand, the new terms appearing in 
H^^ may open up possibilities of new phenomena. In 
the isotropiCj-ajitiferromagnetic case Jx = Jy > 0, Momoi 
and TotsukaEil have shown that correlated hopping pro- 
cesses niib]^,bii" induce the formation of bound states of 
two triplons. In the ferromagnetic case of our main inter- 
est, however, we will argue that the formation of bound 
states is rather unlikely, and that the second-order terms 
in do not change the essential physical properties of 
the system. Instead, the improved accuracy in the es- 
timation of physical quantities allows us to make more 
quantitative comparison with the exact diagonahzation 
results of Sec. 

We first calculate the energy of a single-triplon exci- 
tation. By restricting to the single-triplon sector, we 
can neglect all the terms of the forms, nunni , nfl(l — 
njir)n]^'i and npb^j^ibpjn . The resulting Hamiltonian con- 
tains only the chemical potential /i and the nearest- 
neighbor, second-neighbor, and third-neighbor hoppings, 
whose amplitudes are given respectively by i, —2t', and 
—t'. These parameters are given by 



Jl J2 



2 I J 



8J 
(14) 

The unit cell is again reduced to a single dimer; however, 
this occurs only in a single-triplon sector. This Hamilto- 
nian leads to a triplon dispersion. 



fife = — M + 2f(cos ki + cos fc2) 

- 4t'[cos(fci -I- fc2) -I- cos(fci - ^2)] 
0+' fr^^c ou. I ^^cO^,_^ 



- 2t'(cos2fci + cos 2^2), 

with ki, := k ■ {y = 1,2). Using kx 
ky -.^ k ■ y, this is rewritten as 

Efc — ~IJ,+4:t'+4:tcoskxCosky~16t' cos^ kxCos^ ky. (16) 



(15) 
k ■ X and 
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FIG. 9: (Color online) Contour maps of A/i//ic2 and Ae//ic2 
in the DS-II regions, where the magnetization process shows 
a smooth increase as in Fig. ^(b). Here A/i :— hc2 — hd is the 
width between the starting and ending of the magnetization 
process, and Ae is the width of the triplon dispersion. The 
contours are shown only inside the pink dotted Unes, where 
Eq. ( |C5| ) is vahd. The lines of Ah/hc2 ~ and 1 corre- 
spond to transition hnes to the DS-I region and the stripe 
phases, respectively. The transition lines obtained from ex- 
act diagonahzation analysis of Sec. are also shown (by the 
same symbols as in Fig. ^). Fit with the experimental data 
of (CuCl)LaNb207 yields estimat es o f J^/J and Jy/J, as in- 
dicated by the asterisks (see Sec. VII). 



This dispersion has the minimum energy 

Emin = -/i - 4|t| - 12t', 

and the band width 



Ae = 8|t| 



4|J-|(1 + ^ 



(17) 



(18) 



(if — J+ + 2|J_| < J). The triplon condensation field 
hci in the DS-II regions is determined by the condition 
Emin = 0, leading to 



J 



(19) 



The saturation field hc2 is calculated exactly in Ap- 
pendix Using hc2 = J + 2J+ + 2|J_| in Eq. (|c|), 
which is valid in the coUinear stripe regions of the clas- 
sical model, the width Ah between hd and /ic2 is deter- 
mined as 

Ah :^ h,2 - hd = 2J+ + 4| J_| + ^'^+ ^J^''^' ■ (20) 

In Fig. ^ we display the contour maps of Ah/hc2 and 
Ae/hc2 calculated in this manner. Since Ah = Ae — 
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in the dccoupled-dimer limit Jx = Jy = 0, the two 
parameters indicate to what extent the system is sepa- 
rated from this hmit. However, the two parameters show 
shghtly different behaviors. The contours of Ah/hc2 and 
Ae//ic2 are approximately parallel to the lines Jx = 3Jy 
(or 3Jx ~ Jy) and Jx ~ Jy, respectively. The lines of 
Ah/hc2 — corresponds to the transition between the 
DS-I and DS-II regions. The lines of Ah/hc2 = 1 corre- 
sponds to the quantum phase transitions into magneti- 
cally ordered states associated with the triplon conden- 
sation. These boundaries agree well with the exact diag- 
onalization results plotted together in Fig. ^. This means 
that our perturbative calculation likely gives a good es- 
timation of hci up to rather large values of \Jx\/J and 

\jy\/j. 

Next, we discuss the role of the correlated hopping 
processes nub^^ibun , which appear in the second-order 
effective Hamiltonian. In the isotropic antiferromagnetic 
case Jx = Jy > 0, Momoi and Totsukatil have shown 
that these processes induce the formation of bound states 
of two triplons. Namely, two triplons pair together and 
achieve a lower energy than two independent triplons. In 
this the magnetic field is increased, the bound 

states condense-bofore the triplons do, leading to a bond- 
nematic order.E3E3 It is interesting to look for this pos- 
sibility in the ferromagnetic case of our interest. We 
performed exact diagonalization analysis of the original 
spin model on various points in the ferromagnetic 
case Jx, Jy < 0, and searched for the existence of bound 
states. However, we found no indication of the formation 
of bound states. This rpsjult can be interpreted from the 
study of Schmidt et al. .E3 These authors have analyzed a 
simple hard-core boson model on the square lattice con- 
taining a single-particle hopping t, a correlated hopping t, 
and the nearest-neighbor repulsion V . They have demon- 
strated numerically that the correlated hopping process 
enhances the tendency towards phase separation as well 
as boson pairing; a jump of the magnetization [Fig. ^(c)] 
in the DS-L-iiegion indeed aacurs as a result of the phase 
separation.Ej It was shownEj that the correlated hopping 
process induces the boson paring only when a modest 
magnitude of > is introduced. In the ferromagnetic 
region Jx, Jy < of the model (||), the nearest- neighbor 
interaction is attractive (V < 0), and hence phase sepa- 
ration and a standard superfluid state are more likely to 
occur than the boson pair-condensed state. 



V. EXACT DIAGONALIZATION 

The strong coupling expansion in the previous section 
has provided controlled analyses for small \Jx\/J and 
\Jy\/J. To investigate the ground-state phase diagram 
of the model in a wide parameter space, we have per- 
formed exact (Lanczos) diagonalization calculations in 
finite-size clusters. The shapes of the clusters are shown 
in Fig. where Ng is the total number of spins in the 
system. In every cluster, periodic boundary conditions 
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FIG. 10: Clusters with A^^ = 16, 20, 4 x 6, and 6 x 4 used for 
exact diagonalization analyses. 

are imposed at the edges. The clusters with iV^ = 16 and 
Ns = 20 are invariant under simultaneous operations of 
90° lattice rotation and interchange of Jx and Jy . There- 
fore, these clusters treat Jx and Jy equally, and are less 
biased in studying the ground-state properties than the 
clusters with Ng = 24, which do not have this invariance. 
However, to investigate the size dependence of the cal- 
culated quantities, we analyze the Ng = 24 clusters as 
well. 



A. Magnetization process 

We first investigate the magnetization process in the 
dimer singlet phase. The magnetization process is deter- 
mined by calculating the ground-state energy E{M) as 
a function of the total magnetization M := Sf as in 
Fig. |ll]. To determine the magnetization M for a given 
field h from this figure, we first find the line of slope h/ J 
that touches the curve in such a way as to minimize the 
vertical axis intercept. The magnetization is then found 
at the touching point of the line and the curve. A con- 
cave curve, as in Fig. |ll|(a), leads to a jump of M as in 
Fig. |(a), with 



E{Ng/2) - E{Q) 
Ns/2 



(21) 



A convex curve, as in Fig. IT^fb), leads instead to a 
smooth increase of M as in Fig. E|(b) with 

= E{l)-EiO), hc2 = E{Ng/2)-E{Ns/2-l). (22) 
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FIG. 11: Ground state energy E{M) as a function of the 
magnetization M for (a) Jx = —1 and Jy — —0.5 and (b) 
Jx ~ —1 and Jy = —0.2. A concave curve in (a) leads to 
a jump in the magnetization as in Fig. ^(a), while a convex 
curve in (b) leads to a smooth increase as in Fig. W(b). Broken 



lines indicate the variational estimates using the ansatz (Bl) 



A convenient way to judge whether the curve is con- 
cave or convex is to examine which of he and hd defined 
above is smaller. By plotting he and hd as functions of 
model parameters {Jx/ J or Jy/J) and finding the cross- 
ing point, one can determine the boundary between the 
DS-I and DS-II regions, as indicated by circular symbols 
in Fig. |. Here, the data for = 16, 20, and 6x4 
are shown by symbols of different sizes in order of N^. 
The obtained boundaries show little dependence on the 
system size Ns- In Fig. ^ we also plot the estimate of 
the energy using the Hartree variational ansatz (Bl). We 
find that this ansatz can capture the qualitative shapes 
of the curves. 

As Jx or Jy increases, the spin gap of the system be- 
comes smaller and eventually vanishes, leading to mag- 
netically ordered states. The transition points to the fer- 
romagnetic state can be determined accurately by finding 
the level crossing between the singlet ground state and 
the fully polarized state. These transition points are plot- 
ted by square symbols in Fig. |[ The obtained boundary 
agrees verAi,well with the classical ferromagnetic phase 
boundary.tJ The transitions to the collinear stripe states 
cannot be detected by a level crossing, and we instead 
calculate the fidelity susceptibility, as presented below. 



FIG. 12: (Color online) (a) Fidelity susceptibiUty xf versus 
A := J^ I J for J+ = 0. (b) Extrapolation of the peak position 
of XF into Ns ^ oo using Eq. {Um. 



B. Transition to the stripe phases 

The fidelity susceptibilityii'0 shows a peak at the 
transition point, which we use to detect the phase bound- 
ary. Since the phase boundaries are parallel to the line 
Jx = Jy in the qualitative phase diagram of Fig. ^, it is ef- 
ficient to probe the parameter space in the direction per- 
pendicular to this line. Therefore we fix J+/ J and change 
the parameter A := J-/J, where J± :— (Jx ± Jy)/'^ as 
defined in Eq. (^). The fidelity between two parame- 
ter points A and A' is defined as an overlap between the 
ground states. 



F(A,A') = |(*(A)|*(A'))|. 



(23) 



When the two points are close enough, this quantity has 
an expansion 

F(A,A + <5A) = 1-^XF + ..., (24) 

which allows us to define the fidelity susceptibilityil 

lnF(A,A + ^A) 



Xf(A) = -2 lim, 



(25) 



The result for J_|_ = is presented in Fig. |lj(a). We ob- 
serve a peak in xf(A), which indicates a phase transition; 
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this peak shifts graduahy and grows sharper with increas- 
ing Ng. The peak position Ac(A^s) in ^ finite-size system 
is expected to approach the transition point Ac with in- 
creasing Ng. In Fig. |l2|(b), we extrapolate thCrJata of 
^ciNs) into Ns ^ oo, assuming the scaling form&Zl 



(26) 



with V « 0.711 for an 0(3) quantum phase transition, 
corresponding to triplon condensation. 

The obtained stripe phase transition points are plot- 
ted by diamond symbols in Fig. ^ The phase bound- 
ary crosses the vertical axis at (J^/ J, Jy/ J) ~ (0,0.63). 
When Jx ~ and Jy > 0, the present model is equivalent 
to a Heisenberg antiferromagnet on the honeycomb lat- 
tice with spatially anisotropic exchange coiralings. Quan- 
tum Monte Carlo analysis of such a modeled has given a 
precise estimate Jy/J « 0.576 for the transition point. 
Our result roughly agrees with this, although it slightly 
overestimates the range of the dimer singlet phase. 

The classical phase diagram in Fig. ^ suggests that 
there can be a narrow spiral phase between the ferromag- 
netic and coUinear stripe phases. Since incommensurate 
spin correlations existing in the spiral phase are difficult 
to treat in the small system sizes available to exact di- 
agonalization, we do not analyze the existence of this 
phase within exact diagonalization. Instead, we analyze 
this phase using the Schwinger-boson mean field theory 
in the next section. 



VI. SCHWINGER BOSON MEAN-FIELD 
THEORY 



Here, the local constraint of the form 



(28) 



is imposed so that the spin quantum number S is given by 
the number of bosons per site. The a, /? spin labels run 
over t, i, and cr° with a = x,y, z are the Pauli matrices. 

Heisenberg terms are quartic in the Schwinger bosons, 
so wc apply a mean-field decoupling to obtain a quadratic 
Hamiltonian. The Heisenberg term is given by 



S. -S 



1 



— ^ "l-iiy"pe"ifj."t'^"jp''je 
apvpe 



(29) 



To decouple this, we use the following Pauli matrix prod- 
uct identities to decouple in an attractive channel for each 
possible sign of the interaction: 



(30) 



Here, e^^ is the totally antisymmetric tensor defined with 
Efi = +1. This gives the relations for ferromagnetic and 
antiferromagnetic interactions, 

^S, • S, = -^B.jBI + S{S + 1), (31) 

S.-S, = -^44+^2^ (32) 

where we have defined the operators 

4=E^L^^- (33) 



In this section, we analyze the Heisenberg model in- 
troduced in Eq. (|^) using a Schwinger boson mean-field 
theory.E3c3 The mean-field method can cover a wide 
range of circumstances where our other methods may 
be limited. Unlike the strong-coupling expansion, this 
method is not limited to small Jx/ J and Jy/J- The 
incommensurate spiral state can be easily found by this 
mean-field theory, in contrast with exact diagonalization. 
Furthermore, the theory can determine the splitting of 
the X- and y-spiral state degeneracy due to quantum 
fiuctuations. Finally, one can control the strength of 
quantum fiuctuations by varying the "spin magnitude" 
S, connecting the 5 3> 1 limit of semi-classical ordering 
to the low-spin limit of 5 = 0.5. We consider the case 
h = Q for simplicity. 



A. Theory 

The Schwjager boson representation of the spin on site 
i is given by 



(27) 



A,. 



E ^ia^afibjp. 



(34) 



Wc enforce the boson number constraint by introduc- 
ing a Lagrange multiplier term 



(35) 



We treat this on average by taking A.; constant on each 
of the four sublatticcs. 

Wc perform a mean-field decoupling of the terms given 
by (^). For instance. 



At, A, 



lAijl"^, 



(36) 



with Aj 



(A 



Assuming the translational invari- 



a/3 



ance, there are eight B mean fields, four for the nearest- 
neighbour interactions in the x direction, four for the y 
direction, and two A mean fields for the diagonal anti- 
ferromagnetic interactions. Following this, we perform a 
Fourier transform defined by 

= E e*-'"-+'"^^)&r,;..^. (37) 
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X = A,B,C,D labels a position in a unit cell as in 
Fig. 0(b), so we rewrite the label i — > (r,-'^)- -^uc is 
the number of unit cells in the system. We have used the 
A site position r to label a unit cell, and Srx to represent 
the position of each site relative to the A site, 

5rA — 0, SrB = x, Src — y, Sr^ = x + y. (38) 

At this point, we assume the chemical equivalence of 
the four sites in the unit cell, by taking A; = A the same 
for all sites in the lattice. In conjunction, we use an 
ansatz for the mean fields consistent with the choice of 
chemical potential. Namely, we take the i-direction B 
to be equal (i?x), and do the same for the y-direction B 
(By). The two A mean fields are between sites A and D 
(Aad), and B and C (Abc)- Furthermore, we take all 
mean fields to be real. We take the two A mean fields 
to have equal magnitude, so that A ad = A and Abc = 
±A gives two possible choices for the relative signs of A. 
With these assumptions, the mean field theory can still 



describe all the relevant phases expected to appear in the 
model. 

We are left, in the case of a gapped dispersion, with 
the mean-field Hamiltonian 

nMF ^Hc + J2 bk^kbk, (39) 

k 

where 

Hc = - 2J^Bl - 2JyBl + JA^ - A(8S' + 4), 

bk =(^kAti ^kBtj ^kCti ^kDt; 

^-kA^' ^-kSi' ^-kCi' ^-k_Di 

). (40) 

The matrix Hk is given by 



I A JxBxCoskx Jy By cos ky 

JxBx cos kx A 

A JxBxCoskx 

Jy^v *^os ky JxBx cos kx A 



JyBy COS ky 



J y -By COS ky 



(42) 



Ua 



Aoe 






iABce'^"^ 





-^Abcs^'^''''^''' 





-lA 



(43) 



with kx = k ■ X, ky = k ■ y, and Aadt Abc as defined 
above. 

The quadratic Hamiltonian in ( |39| ) is diagonalized by 
the Bogoliubov transformation bk = ^fc7ki with 7 de- 
fined in the same manner as b: 

7k = (7kAt, 7kBt, 7kCt, 7kZ?t, 

"^-kAJ,' ^IkSJ.' ^IkCJ,' '>'-k_D^) ■ 

(44) 

This is a canonical transformation, where the operators 
in 7 preserve the bosonic commutation relations of the 
operators in b. These commutation relations are given 

by 

[7k,7k'] ='5kk'??, (45) 
where the 8x8 matrix ?/ can be written in 4 x 4 blocks: 

„^(J_°,), (4.) 

The proper commutation relation are obtained by tak- 
ing the columns of Zk to be the eigenvectors y^k (with 



eigenvalues Wmk) of ry'Hk- Here, m = 1, . . . , 8 labels the 
eight eigenvectors. 

When the Hamiltonian Hk is positive definite, the cor- 
responding spinon dispersions IwmkI are gapped. This 
leads to a disordered ground state. In the case of a gap- 
less dispersion, ryHkoym'ko = for some to', so that both 
"Hko and ??7^ko have a zero eigenvalue at the dispersion 
minimum kg. Condensation of such zero-energy bosons 
gives rise to a magnetically ordered state. To describe the 
condensation, we replace the operators bko with macro- 
scopic constant values, 

^™'ko = X! '^"'•'o V^y™'ko 1 (47) 

m' m' 

where yti/koym'ko — 1- With the condensate contribu- 
tion, the diagonalized Hamiltonian is 

8 

771K I / mk 

k ni—1 

+ EE^»'ko^k„x™.ko. (48) 

ko rn' 
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Spin ordering is found at twice the spinon minimum 
wavevector, 2ko, governing relative spin orientation be- 
tween unit cells. 

Given the diagonalized Hamiltonian, we must solve 
for the mean- field values of A,Bx,By and A, as well as 
any condensate vectors Xm'ko ^-^d associated minimum 
wavevectors kp. It turns out that the solution depends 
only on the total condensate density ^j^^ J2m' |cm'koP- 
We have the mean-field equations 



MFI 



MFI 



MFI 



dA 



OB, 



OB, 



MF 



d{H 



MF 



dX 



0, 



(49) 



which we solve self-consistently for a given set of param- 
eters Jx/J, Jy/J, and S, and the associated values of m' 
and kp. 




1 



Jx/J 



B. Results 

In the semi-classical limit S oo, we recover the clas- 
sical phase diagram shown earlier in Fig. ^. We present 
the phase diagram as a function of Jx/J and Jy/J for the 
cases of S' = 0.5 (Fig. 1) and fi- = 0.15 (Fig. |l|). Since a 
mean-field theory is expected to underestimate quantum 
fiuctuations, it can be instructive to look at spin values 
smaller than the actual case. The Schwinger boson mean- 
field theory finds several magnetically ordered phase, and 
one disordered phase. The magnetically ordered phases 
are the ferromagnetic, spiral and stripe phases seen in 
the classical limit, while the disordered phase features 
isolated dimers. Below we summarize how each phase is 
described in the Schwinger boson formalism. 

Ferromagnetic Phase. In this state, the antiferromag- 
netic order parameter A = 0, while both ferromagnetic 
order parameters take their maximum values, Bx = By = 
2S. The spinon minimum wavevector kp =0, as is the 
spin ordering wavevector, so that the spins are fully po- 
larized throughout the entire lattice. 

Stripe Phase. In this state, the ferromagnetic order 
parameter associated with the smaller of Jx and Jy is 
zero. The ordering wavevector is 0, and a stripe magne- 
tization pattern is found in the direction of the larger fer- 
romagnetic interaction. Depending on the relative signs 
of the A parameters, the spinon minimum wavevector is 
either (0,7r) or 0. for Jx > Jy, both leading to the same 
ordered state of same energy. 

Spiral Phase. This is an incommensurate magneti- 
cally ordered state with a nonzero ordering wavevector 
in the direction of smallest ferromagnetic interaction: 
ko = {kx , 0) for Jy > Jx- All of the mean-field parame- 
ters are nonzero. There are two inequivalent spinon min- 
imum wavevectors ±ko found in the spiral state. They 
both lead to the same type of spiral order, described in 
Sec. III. The relative sign of the two A mean fields yields 
either the x-spiral or y-spiral ordering. We find that 



FIG. 13: Schwinger boson phase diagram for S — 0.15. Solid 
lines indicate second-order transitions, while dashed lines in- 
dicate first-order transitions. The dotted line separates the 
two spiral phases. The disordered dimer state (DS) expands 
significantly compared to the S = 0.5 case in Fig. ^. The 
stripe state expands moderately, pushing the ferromagnetic 
(FM) state to moderately larger Jj, and Jy. In between these 
phases, the spiral state is found in a significantly reduced 
area. With an even smaller S = 0.1, the spiral state will have 
disappeared entirely. 



the classical degeneracy between x- and y-spiral order- 
ing is broken by quantum fiuctuations, determining the 
ordering direction. For Jx > Jy, the y-spiral state has 
a lower energy, while the a;-spiral state has lower energy 

for Jy > Jx- 

Disordered Phase. The disordered phase is gapped, 
with no boson condensate, and has no ferromagnetic cor- 
relations: Bx = By = 0. The antiferromagnetic dimer 
order parameter A ^ Q, leaving a decoupled dimer state. 
This is a mean field description corresponding to the 
dimer singlet state. 

We consider the effect on the phase diagram of low- 
ering spin from the large- S" semiclassical limit, by com- 
paring the classical, S = 0.5, and S = 0.15 phase di- 
agrams shown in Figs. |[ |^ and ^ respectively. As 
S decreases, the disordered phase appears and expands 
around Jx ^ Jy = 0, pushing out the magnetically or- 
dered phases to larger Jx and Jy. Of particular notes is 
the shift of the ferromagnetic phase boundary to larger 
Jx and Jy in contradiction with the exact boundary for 
Jx = Jy = J seen in Appendix ^ Similarly, the stripe 
state boundary shifts to larger Jy for Jx > Jy, and larger 
Jx for Jx < Jy. Between these aforementioned phase 
boundaries, we find the spiral state, which consequently 
gets pushed out to larger Jx and Jy. This state shrinks 
as S decreases, having almost disappeared at = 0.15, 
as seen in Fig. O. In contrast with the case of the purely 
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antiferromagnetic Shastry-Sutherland lattice^ we see 
no other disordered states than the dimer one down to 
5 = 0.1. In particular, there are no short-range ordered 
analogues of the ferromagnetic, stripe, or spiral ordered 
states. 



VII. COMPARISON WITH (CuCOLaNbjOr 

In this section we compare our theoretical results with 
the experimental data on (CuCl)LaNb207. The mag- 
netization process of the compound shows a smooth in- 
crease between Bd = 10.3 T and Bc2 = 30.1 T,a as in 
Fig. ||(b). Therefore, the compound should be located in- 
side the DS-II regions in Fig. |^ within the J-Jx-Jy model. 
Assuming the Lande factor g = 2, the above magnetic 
fields are translated into energy units as 



hci = gfJ'BBci = 1.2 mcV, 
hc2 = gfJ-BBc2 = 3.5 meV. 



(50) 



According jto a recent inelastic neutron scattering 
expcrimentjij the triplet excitations at _B = range from 
approximately 1.2 to 3.0 mcV. Good agreement of the 
lower bound of the excitations with hd suggests that 
the emergence of the magnetization is associated with 
the condensation of single-triplon excitations. Indeed, a 
Bose-Einstein condensation of some sort of bosonic ex- 
citations haV|e been observed in the field range between 
Bel and Bc2-^ Here we do not address the possibility of a 
condensation of bound magnons proposed in Refs. ^|J3^, 
since it does not seem to occur within the J-Jx-Jy model 



(except for the case J^, > 0) as discussed in Sec. |IVC . 

Below we first discuss the experimental data within 
the J-Jx-Jy m odel, using the second-order perturbative 
results of Sec. IV C . Then we discuss possible effects of 
other exchange couplings J2b and J'^ shown in Fig. 0(a), 
within the first-order perturbation theory. 



J-Jx 



Using the above data of the magnetization process and 
the triplet excitations, the two ratios calculated in Fig. ^ 
are given by 



A/i//ic2 = 0.66, Ae//ic2 = 0.52. 



(51) 



These can be used to fix the values of ( J^,/ J, Jy/ J), as 
shown by asterisks in Fig. Adjusting J to give hc2 in 
Eq. (^), we obtain the estimates 

J = 2.2 mcV, JxlJ = -0.08, Jy/J = 0.29. (52) 

Another parameter set where the values of Jx and Jy are 
interchanged is also possible. In these naive estimates 
within the J-Jx-Jy model, one of Jx and Jy is ferromag- 
netic while the other is antiferromagnetic. This is in con- 
trast to both ferromagnetic values Jxr-Jy < predicted 
in the electronic structure calculation.E3 




(7tf2,7C/2) 



(n/2,%/2) 



FIG. 14: (Color online) Triplon dispersions ( p4| ) in the J- 
Jx-Jy model in the cases of (a) ferromagnetic Jx < and 
antiferromagnetic Jy > in Eq. ([5^) and (b) both ferromag- 



netic Jx, Jy < in Eq. (| 



Black solid and broken lines 



correspond to ek,- and ek.+, respectively. These are plotted 
along a triangle connecting (7r/2, 7r/2), (0,0), and (7r/2,0) in 
the Brillouin zone. Blue horizontal dotted lines indicate the 
energy range of triplet excitations (from 0.34/ic2 to 0.86/ic2) 
observed in Ref. PL3. To compare with the magnetic layer of 
(CuCl)LaNb207 in Fig. |l](a), we have replaced k — > ka, where 
a ~ 3.9A is the distance between neighboring Cu^^ ions. 



We note, however, that the estimation depends on 
the interpretation of the experimental data. In addi- 
tion to the excitations ranging between 1.2 and 3.0 meV 
above, an earlier inelastic neutron scattering cxpcrimcntti 
also detected excitations centered around 5.0 meV (w 
1.43ft.c2) with small scattering intensities. In Rcf. |^, these 
excitations are interpreted as bound states of triplets. If 
we regard that these also originate from triplet excita- 
tions, the range Ae is much larger. Then in Fig. |[ the es- 
timates of {Jx/ J, Jy/J) can change, keeping the relation 
Ah/hc2 = 0.66, to the case of ferromagnetic Jx, Jy < 0. 
For example, in the parameter set 

J = 4.4mcV, J^/J= -0.82, J^/J = -0.10, (53) 

the triplon excitation energy in Eq. (^6|) ranges up to 
l.32hc2- It is expected that by increasing \Jx\ further, 
the energy upper bound would reach 1.43/ic2 observed in 
experiment. We do not discuss a larger- jJrl region, since 
our perturbative result would be less reliable there. 

To further discuss the consistency between the present 
model and (CuCl)LaNb207, it would be useful to look 
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(71/2,71/2) 



(0,0) 



(7C/2,0) (7C/2,7t/2) 
ka 



(0,7C/2) 



(0,0) 



FIG. 15: (Color online) Triplon dispersions ( pq ) in the J-Jx- 
Jy-J2b-JA model. Compared to Fig. ^ the path (7r/2, 7r/2) — ^ 
(0,7r/2) (0,0) is added, since the triplon dispersions are 
anisotropic in the and ky directions in the presence of 
nonzero J26 or J4. 



at the detailed shapes of the triplon dispersions. Here 
we plot the triplon dispersions using the perturbative re- 
sult in Eq. ([iq). When more detailed information of the 
triplet excitations is provided from experiments, these 
plots can be used to test the J-Jx-Jy model and to de- 
termine the signs of Jx and Jy. Equation (^6|) con- 
tained only a single band, since the single-triplon hop- 
ping problem had the same periodicity as the square lat- 
tice (of dimer centers). To discuss the physical excita- 
tions, however, one needs to take into account the fact 
that there are two dimers in a unit cell of the original 
lattice, and to fold the Brillouin zone in such a way as 
to identify k, k -|- (tt, 0), and k -I- (0,7r). Consequently, 
the first Brillouin zone is given by a square ranging over 
— 7r/2 < kx^y < tt/2. After the folding, there are two 
dispersions expressed as 



16t' cos^ kx cos^ k 



ek ± ~ + 4:t' ± 4:\t\C0S kx cos ky J^^L, ^^.^ n,^; ^.WO n,y. 

rrn ^^^^ 

These dispersions are plotted in Fig. |j for the pa- 
rameter sets in Eqs. ^ and (||). While Fig. |l|(a) 
shows a nearly symmetric shape centered around 0.66/ic2, 
Fig. |lj(b) shows a highly asymmetric shape. 



B. 



J-Jx-Jy-J2b-J4 model 



In addition to J, Jx, and Jy, the model of Ref. |IJ also 
contains some other exchange couplings such as J2b and 
J4, shown in Fig. 0(a). Here we discuss their effects on 
the triplon dispersions. 

In the first-order perturbation theory, the additional 
couplings above introduce triplon hopping and interac- 
tion terms between the dimers they connect, similar to 
Eq. (^. Using the effective Hamiltonian, triplon disper- 
sions at ft, = are calculated as 



(55) 



with 



JL 



cos4fcx cos2ky, 



(56a) 



(5k = —{Jx — Jy) COS kx COS ky + J2b COS 3kx cos ky, 

(56b) 



J' 

i?k = sin 4^2; sin2fcy. 
Around k = 0, these dispersions are expanded as 



(56c) 



ek,± -J±|AJ|-:|t 



lAJl 



4(sgn AJ)J26T44 



(57) 



with A J := — Jx + Jy + J26. This dispersion implies that 
the minimum of the triplon excitation energy is given by 



(58) 



.o.- = J-|AJ|-|, 
if the following inequalities are satisfied: 

^ + 4(sgn AJ) J2fa -t- 44 > 0, ^ + 4 > 0. 

(59) 

Eauation (M) can be used to determine the critical fields 



(60) 
(61) 



J' 

hc2 = J+Jx + Jy + J2b + Y + l^-^l- 



The electronic structure calculation of Ref. [l^ pro- 
duced the estimates 



Jx/J = -0.39, Jy/J = -0.38, 
J2b/J = -0.14, 4/J = 0.18. 



(62) 



These parameter values satisfy the condition (|5^) and 
give hci/J = 0.78 and hc2/J = 0.31. We find hd > hc2, 
which means that the system is in fact inside the DS- 
I region, and the magnetization jumps as in Fig. ^a). 
We confirmed the occurrence of a jump by directly per- 
forming an exact diagonalization for the J-Jx-Jy-J2b-Ji 
model with the values of Eq. (^2|). In order to be con- 
sistent with the smooth, increase of the magnetization 
observed in experiment the estimates in Eq. (^2|) must 
be modified. 

Given the difficulty in uniquely determining the ex- 
change couplings on the basis of the current experimen- 
tal data, we here discuss possible effects of J26 and J4 
on the triplon dispersions, which could be tested in fu- 
ture experiments. We add small values of J2b and J4 to 
the parameter set in Eq. (p3[) , and plot the triplon dis- 
persions Eq. (H) in Fig. (^. Compared to Fig. |l|(b), 
the dispersions show oscillating behaviors coming from 
the long-range triplon hopping induced by J2b and J4. 
These oscillating behaviors can be used as fingerprints 
for the existence of these couplings. 
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VIII. CONCLUSIONS 

Motivated by recent experiments on (CuCl)LaNb207, 
we have studied a spin-^ model (0) of ferromagneti- 
cally coupled dimers on the distorted Shastry-Sutherland 
lattice. Using the three different approaches (the 
strong-coupling expansion, exact diagonalization, and 
Schwinger boson mean field theory) , we have determined 
the ground state phase diagram of this model, as shown 
in Figs. ^ and ^. We have shown that in the dimcr sin- 
glet phase appearing for |Ja;|,|Jy| < J, the magnetiza- 
tion process depends crucially on the spatial anisotropy 
of Jx and Jy. In the DS-I region with weak anisotropy, 
the magnetization shows a jump, as in Fig. |(a). In the 
DS-II regions with strong anisotropy, the magnetization 
smoothly increases after the spin gap closes at a certain 
magnetic field, as in Fig. 11(b). When or |J^| > J, 
quantum phase transitions to various magnetically or- 
dered phases (ferromagnetic, collinear stripe, and spiral) 
occur. These magnetic phases also appear in the classi- 
cal limit of the model as shown in Fig. ^. Providing a 
notable difference from the classical case, the Schwinger 
boson analysis has demonstrated that quantum fluctua- 
tions split the classical degeneracy of two kinds of spirals 
(x- and j/-spirals). 

We have compared our theoretical results with the ex- 
isting experimental data on (CuCl)LaNb207. A smooth 
magnetization process observed in the compoundcl sug- 
gests that, at least within the J-Jx-Jy model, the com- 
pound should be located in the DS-II region, and that 
there should be substantial anisotropy in and Jy. 
Comparing the second-order perturbative results with, 
the experimental data for ±lie magnetization processci 
and the triplon excitations J13 we found that one of Jx 
and Jy may be ferromagnetic while the other may be 
antifcrromagnctic. This is in contrast to both the ferro- 
magnetic estimates in the previous electronic structure 
calculation. Ej However, we leave open the possibility of 
both ferromagnetic Jx and Jy, given the possible exis- 
tence of higher-energy triplet excitations. We plotted 
the triplon dispersions for different signs of Jx and Jy. 
These plots can be used to test the J-Jx-Jy model fur- 
ther and to determine the signs of Jx and Jy, when more 
detailed information of the triplon excitations is provided 
from experiments. We have also discussed the effects of 
other exchange couplings J26 and J4 [shown in Fig. 0(a)] 
on the triplon dispersions. 

Our analysis of the distorted Shastry-Sutherland 
model may provide a useful starting point for under- 
standing the wide range of physics in the layered copper 
oxyhalide family. A quantum phase transition from a sin- 
glet ground state to a collinear stripe state observed by 
chemical substitutioicHiij may be interpreted as the tran- 
sition from the DS-II region with Jx < Jy to the stripe 
state with a propagation vector q = (0, tt) in Fig. ^. Here 
we have excluded the case of Jx > Jy, since the stripe 
state with q = (tt, 0) appearing in this region is in fact a 
Neel state with q = (tt, tt) in the original square lattice 



of the compound in Fig. |T|(a). The ion-exchange method 
used to synthesize the layered copper oxyhalide family 
has a remarkable flexibility of changing the constituent 
atoms in various ways. We expect that the jump of the 
magnetization in the DS-I region and a spiral order with 
a unique structure as in Fig. ^ will also be observed in 
future experiments on related compounds. 
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Appendix A: Exact ground states in the isotropic 

CclSG rJx — Jy 



Here we discuss the exact ground states of the Hamil- 
tonian (|^) in the isotropic case Jx = Jy{= J')- The 
discussion of this sccticp, is based on the argument of 
Shastry and Sutherlandllfl and its extension to the case 
of ferromagnetic J'. We decompose the Hamiltonian H 
into triangular parts and solve a 3-spin problem as in 
Fig. @(a). The 3-spin Hamiltonian is given by 



i?A = J'{Si-S2 + S2-S^) + ^Si-S^~^{Sl + SI). (Al) 



Note that here the Zeeman coupling with the magnetic 
field is introduced only for the first and third spins and 
not for the second spin. The Zeeman term for the second 
spin is included in the Hamiltonians of different triangles. 
In general, the decomposition of H to triangular pa rts is 
not unique. However, definition of H/^ usfid in Eq. ( Al), 
introduced by Shastry and Sutherland,llJ will be useful 
in the following discussion. 

When /i = 0, Hi^ is S'J7(2)-symmetric, and its cigcn- 
states consist of a quadruplet \q^) (with /i = ±i,±| 
labeling S^^^) and two doublets. Because H^^ is also sym- 
metric under the permutation of the sites 1 and 3, the 
doublets are classified into ones symmetric and antisym- 
metric with respect to this operation, \d+i) and \d'i) 

respectively. When h ^ 0, the SU{2) symmetry of Ha 
is reduced to U{1). Then \q±^) and |c?-|-i) are mixed to 

form new eigenstates, \q±i) and \d^i). The eigenstates 
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arc summarized as follows: 



i9+i> = itn), 

|g+i(7+)> = 



(A2a) 



1 



2 + 
1 



:(im> + i;n> + 7+int)), 



27i 



(A2b) 

(7+im)+7+i;n)-2int», 



V2 



(im>-i;n», 



(A2c) 
(A2d) 



and the other eigenstates with /i = —1/2 or —3/2 arc 
obtained by reversing all the spins and replacing 74- with 
7_ . The parameters 7± arc given by 



The corresponding eigenenergies are given by 



8 2' 



8 " 7±^ 
3J 



(A3) 

(A4a) 
(A4b) 
(A4c) 
(A4d) 



Crucially, the two states \d'i) have the same eigenen- 

2 

ergies and are written as product states of a singlet on 

the sites 1 and 3, and an isolated spin state on the site 

2. The exact dimer singlet state I^'ds) in Eq. (||) can be 

written in terms of \d'i) on a given triangle. The state 
* 2 

I^Ds) is therefore an eigenstate of the triangular Hamil- 

tonian 7?a- If are the ground states of the tri- 

2 

angular Hamiltonian, the dimer singlet state I^'ds) will 

also minimize the energy of the triangular Hamiltonian. 

Consequently, I^'ds) becomes an exact ground state of 

H, since it minimizes the local energy of every triangular 

Hamiltonian. We note that the degeneracy of \d',i) is 

2 

not split by a magnetic field h, which comes from the 

specific form of the Zeeman coupling in Eq. ( |AlD . 

Assuming J > and /i > 0, the ground state of i?A is 

summarized in Fig. 16. For ft, = 0, \d'i) are the ground 
' — ' * 2 

states of i?A when 



1 J' 1 

^<7<2- 



(A5) 



These states remain the ground states of H/^ up to a 
certain magnetic field, indicated by solid lines in Fig. |l^. 
In this regime, the dimer singlet state I^'ds) is an exact 
ground state of H with eigenenergy — |A^s J. 
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FIG. 16: The ground state of the 3-spin Hamiltonian Ha- 
The solid lines indicate the boundaries on which the 3-spin 
ground state changes. The exact ground states of H derived 
from the 3-spin solution are indicated in the parentheses. For 
J' > 0, the exact dimer singlet state |^'ds) continue to be the 
ground state of H beyond the solid line (see the text) . To indi- 
cate the range where |^'ds) actually remains the ground state, 
we show broken and dotted line. These lines respectively in- 
dicate the magnetic fields on which two- and one-triplon ex- 
citations condense according to the perturbative calculation 
of Ref. |l|. 



For J' < 0, the ground state of i/A changes to |<?-|-|) 
as we pass the solid line in Fig. |l^. In this regime, the 
ground state of H is given by a fully polarized ferro- 
magnetic state since it contains the state |(7_|_3) locally 
on every triangle. Combined with the above result, the 
present exact argument can complete the phase diagram 
of H in the case of ferromagnetic J'. Specifically, the 
phase boundary is given by /i = J + J'. 

For J' > 0, the ground state of i/A is replaced by 
M+i) beyond the solid line in Fig. |l^. The solid line is 
give by 



A{J + J'){J - 2J') 
4J- J' 



(A6) 



Outside of this region, we cannot use the 3-spin solution 
to predict the ground state of H, since l^+i ) is not fac- 
torizable. However, one can show that I^'ds) is always 
an exact eigenstate of H, regardless of t he values of J, 
J' and h. Therefore, beyond the line (A6), I^'ds) should 
still remain the exact ground state of H until it encoun- 
ters a level crossing with another eigenstate. For h — 0, 
this crossing has been found to occur at J'/ J sji_.0.69 
in an exact diagonalization study of finite clusters.113 For 
/i > 0, the level crossing is associated with a condensation 
of triplon excitations. As shown in Ref. a bound state 
of two triplons condense before one-triplon state does. In 
Fig. |l^, we draw the two lines, derived from the strong 
coupling expansion to third order in Ref. where one- 
and two-triplon condensations occur. 
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Appendix B: Hartree variational state 

Here we discuss an alternative derivation of the p1i « 
diagram in Fig. |[ using the Hartrcc variational stateOo 
The discussion goes in parallel with Sec. IV, and some 
notations are common with this section. The Hartree 
variational state is given by a product of local states on 
dimers: 



R 



^> = n ( cos + e^*« sin 



(Bl) 



This state smoothly connects between the dimer singlet 
state {Or = tt) at /i = and the fully polarized state 
(Or = 0) at a high field. For < Or < ir, this state 
describes a superfluid state of triplons, which breaks the 
U{1) symmetry. We minimize the energy (\E'|i/|4') with 
respect to the variational parameters {0r; 4>r} to find a 
variational ground state. Using the obtained values of 
{^r; (t>R}i one can determine the magnetic structure via 
the relations 



Rl 



'Rl' 



R2 + 



1 



R2/ 



2^/2 



sin 9r, 



Rl/ 



R2l 



^{I+cosBr). 



(B2a) 
(B2b) 



Since I^E") G Fq, we find 

= {^PHP\^) = Eo + {'<if\Hf\'<if), (B3) 

wher e P is the projection operator onto Vq. Since 
can also be viewed as a coherent state of sr 
'i'\Hi^\'9) is precisely equivalent to 



Eq. (§1 
spins, minimizing 



finding the classical ground state of tlic-X2 

cmEacJt 



model Hf. 

The solution to the classical problemlfS'EJo is essentially 
similar to the quantum solution (except for the detailed 
shape of the magnetization process), and leads to the 
same phase diagram of the original model H as in Fig. ^. 
In the entire parameter space, the obtained variational 
ground state has the same periodicity as the original 
model; the values of Or and 0^ arc uniform in each sub- 
set of dimers {T>i or )• Therefore, the ground states 
are characterized by four parameters 0i 2 and 4'i,2 defined 
for I?i.2- In the DS-II region with Jx < Jy, the classical 
ground state of Hf^ is given by a canted antifcrromag- 
netic state with 0i = (j)2+Tr and 9i =62- Using Eq. (B2), 
this leads, in the original model, to a canted stripe state 
with a propagation vector (0, tt) in the transverse com- 
ponent. Similarly, in the DS-II region with > Jy, the 
classical solution of is given by (f>i — 02 and di = 6*2 , 
leading to a canted stripe state with a propagation vec- 
tor (tt, O). The magnetization curve in the DS-II region is 
determined as in Fig. ^(b) , but the curves are all straight 
in the current variational approach. As we increase \Jx\ 
or \Jy\, hci goes to zero, signaling a phase transition to a 
collinear stripe phase. We note that even after the tran- 
sition, the variational ground state at h = remains the 



dimer singlet state I^'ds) with Oi = 82 = tt. This is an ar- 
tifact of the ansatz. In the DS-III region, the variational 
solution for M = A/.s/2 is given by (0i,02) = (0,7r) and 
(tt, 0), indicating th existence of the plateau. 



Appendix C: Single magnon from the polarized state 

Here we consider a single-magnon excitation from the 
fully polarized state |FM) = | ftt ■ • ■) and determine 
the saturation field hc2- We introduce the Fourier- 
transformed basis 



|X,fe) 



1 



(CI) 



where X = A,B,C,D labels a position in a unit cell as 
in Fig. 0(b). We have used the A site position r to label 
a unit cell, and Srx represents the position of each site 
relative to the A site, as in Eq. (|38|). Using the basis 
(|Cl|), we define a 4 x 4 matrix A^(fey by 



Mxx' (fc) = {X, k\H\X\ k) - 5xx' {FM\H\FM) . (C2) 

Here, the energy of the fully polarized state is subtracted. 
To represent the matrix A^(fe) in a compact form, we 
identify X ^ A, B,C, D with fictitious two-spin-i states 
tt, tij iti ii and introduce two sets of Pauli matrices, cri 
and (72. The obtained expression is 



M{k)=J{k)+ [h-^-Jx-Jy 



(C3) 



where 



J'(fc) — Jx(T2 cos kj. + Jyaf cos ky 

+ -^{cTi cos kx + trf sin k^ ) (o-f cos ky + o-| sin ky ) . 

(C4) 

with kx ~ k-x and ky = k-y. By diagonalizing J^{k) and 
obtaining the eigenvalues X^{k) [v = 1,2,3,4), one ob- 
tains four magnon dispersions \u{k) :— X^{k) + h — J/2 — 
uming that the leading instability from the 
fully polarized state is a condensation of single-magnon 
excitations, the saturation field hc2 is found when the 
minimum of the lowest dispersion Ai(fc) touches zero. 
This leads to an expression hc2 ~ ^ + Jx + Jy + Amin, 
where Xnun is the minimum of Ai(fc) over k. 

We now determine Amin- One cap-show, through the 
formulation of Luttinger and Tisza,c3 that this problem 
is precisely equivalent to the classical ground-state prob- 
lem. Therefore, the analysis can be done in parallel with 
Sec. HI. In the collinear stripe phases of Fig. ^, the min- 
imum of the dispersion is found at A; = (0,0). At this 
point, J'(O) is expressed using only mutually commuting 
matrices, af and erf, and thus is diagonalized easily. The 
saturation field is 



J^Jx^Jy^lJx Jy 



(C5) 
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In the classically spiral regime, the minima of the disper- 
sion are fomrd at fc = {Qx,0) and (0, Qy), where Qx and 
Qy are given in Eq. (||) . When one of and ky is set to 
zero, the diagonalization of J^{k) is simple, because one 
of erf and erf becomes a good quantum number leaving a 
single-spin problem in a magnetic field. The saturation 
field is calculated as 

jj^jy /2 1 1 y 

We comment that in the isotropic case Jx = Jy, the 



minimum of the dispersion is obtained along a circle in k 
space, leading to an interesting pseudo-one-dimensional 
density of states of magnons.Eil 

In closing, we stress that our argument for determin- 
ing hc2 in this section is based on the assumption that 
the leading instability from the fully polarized state is a 
condensation of single-magnon excitations. The present 
argument therefore does not apply to the DS-I region, 
where a direct transition occurs from the fully polarized 
state to a singlet state as one lowers the magnetic field. 
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